Dynamics of a trapped domain wall in a current perpendicular to the plane spin valve 
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Abstract 

A study of transverse tail-to-tail magnetic domain walls (DW) in novel current perpendicular to the plane (CPP) spin valves 
(SV) of various dimensions is presented. For films with dimensions larger than the DW width, we find that DW motion can 
give rise to a substantial low frequency noise. For dimensions comparable to the DW width, we show that the DW can be 
controlled by an external field or by a spin momentum torque as opposed to the case of CPP-SV with uniform magnetization. 
It is shown that in a single domain biased CPP-SV, the spin torque can give rise to 1/f-type noise. The dipolar field, the spin 
torque and the Oersted field are all accounted for in this work. Our proposed SV requires low current densities to move DW 
and can simulate devices for logical operation or magnetic sensing without having to switch the magnetization in the free layer. 



PACS numbers: xxx 



I. INTRODUCTION 



The study of domain walls (DW) magnetic (static and 
dynamic )i*2i 3 j4 and transpor^Si^ properties have at- 
tracted much attention recently due to their relevance 
for magneto-electronic nano-device applications. In the 
static case, McMichael and Donahue have shown that in 
magnetic stripes head-to-head DW can be in vortex or 
transverse shaped Klaui et al. have observed these DW 
structures in magnetic thin films and rings 10 . In this 
work, we study the dynamics of DW with magnetization 
mostly restricted to the plane of a nanometer-size thin 
film. DW motion can be a source of low frequency noise 
in magnetic systems and thus hinder any potential appli- 
cations. Their high frequency behavior can however be 
tuned by choosing suitable boundary conditions or ge- 
ometries. Constrained DW oscillations can be high in 
frequency and hence there is a potential for their use in 
nano-electronics as resonators^. 

At least two modes of oscillations have been studied 
in the literature; a Doring-type oscillation and a Winter 
mode oscillationAiiiSiiiiii. The Doring mode is associ- 
ated with translations of the center of 'mass' of the DW 
in an infinite system, while the Winter modes are non- 
zero energy modes that, in addition to rigid translation, 
correspond to propagations along the DW. These modes 
are found by solving the time-dependent Landau-Lifshitz 
(LL) equations*i£J£. 

In the following, we focus only on dynamical properties 
of transverse tail-to-tail DW trapped in stripes with di- 
mensions comparable to their width. This case is proved 
to be the most interesting due to the unusual magneti- 
zation dynamics as opposed to the conventional case of 
DW with widths being much smaller than the film size. 
In particular, we find that the low frequency excitations 
can be reduced in comparison with the usual CPP config- 
uration having a uniform magnetization. Such uniform 
magnetization is susceptible to large fluctuations due to 
spin momentum transfer and gives rise to appreciable 
1 / /-like noise. 17 We investigate DW formed by pinning 



the magnetization in the opposite directions at the edges 
along the easy axis. We show that DW motion can be 
controlled by a current perpendicular to the plane of the 
DW magnetization in contrast with the currently actively 
studied case of DW motion in nano-wires^ 



The LL equation is the basis for the present study. 
Since, we are looking at thermal and current effects, a 
random field and a spin torque term are also added to the 
LL equatior»iS*i2iSfi. We show that in the spin valve (SV) 
geometry suggested here the spin torque can provide the 
force needed to move the wall in a controlled fashion with 
100 times smaller current values than those needed to 
switch the uniform magnetization in a CPP SV. 



The paper is organized as follows. In section [n] we 
describe both the theoretical background and the com- 
putational model needed for our study. In section IIIII 
we discuss the DW structure and we also calculate the 
lowest eigenmodes of the DW using a simple one dimen- 
sional model and compare it to the numerical solution. 
It is shown that if the center of the DW is allowed to drift 
along the easy axis, the contribution to the 1 /f-type noise 
increases. In section ITVl we study the effect of an exter- 
nal field and a CPP current on the SV with and without 
DW. We find that SV with a uniform magnetization can 
be susceptible to unwanted behavior due the spin torque 
driven instabilities that are absent in a SV with a con- 
strained DW. We consider a new CPP geometry with 
a constrained DW between two fixed layers one pinned 
along and the other pinned perpendicular to the direction 
of the current flow. We find that in this CPP structure 
the DW motion can be well controlled with current densi- 
ties which do not lead to the magnetization instabilities. 
It is shown that accurate quantification of demagnetizing 
fields is not essential for a qualitative understanding of 
the influence of the spin torque on the DW. In section 
Ivl we summarize our results. 
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II. THEORETICAL BACKGROUND AND COM- 
PUTATIONAL MODEL 

For a magnetization M with magnitude M s , the LL 
equation for m = M/Af s with a damping in the Gilbert 
form and time normalized by "fM s , where 7 is the gyro- 
magnetic ratio, is given byi& 



dm 



= m x 



,eff 



1 1 \ ^ m 

K(t) - a— 



(1) 



where the effective field h cff includes the exchange inter- 
action, the anisotropy field along the x axis, the demag- 
netization field, the Oersted field, and the spin torque 
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2A. 



V 2 m 



2K 



x(m • x) + hd + plm x m p , (2) 



The damping term is taken to be a = 0.02 in the absence 
of currents and is increased to a — 0.08 in the presence 
of spin torques to account for spin accumulation at the 
normal-ferromagnetic interface. The exchange constant 
A = 1.6 x 10 -6 erg/cm in this study. The random field 
h r (t) is taken to be uniform and Gaussian white at tem- 
perature T, < hr ti (t)hr,j(t') >= 2akT/(jM^V)S i:j 5(t - 
t')i& In the presence of spin torques, the white noise 
assumption — is strictly valid only for frequencies around 



the resonant frequency as shown in Rcf. |20J since we are 
only interested in currents below the critical current, the 
white noise assumption will not alter the qualitative con- 
clusions of this work. The size of the discretized cell is 
taken 2x2 nm 2 in the plane of the film. The inclusion 
of the demagnetizing field is important in DW motion 
studies and hence a numerical treatment is often needed 
to get a quantitative understanding of the dynamics of a 
DW 22 i 23 i 24 . The last term in Eq. His the contribution of 
a spin torque from the pinned layer (PL) m^. The pref- 
actor p is dependent of geometrical parameters and / is 
the current flowing perpendicular to the magnetic multi- 
layers. The p prefactor is dependent on the thickness d, 
the cross section A of the layer and the polarization of 
the current. Assuming perfect polarization of the con- 
duction electrons, with charge e, by the reference layer 
and neglecting the angular dependence dependence, the 
spin torque coefficient is given by 

P= \e\dANir (3) 

In the following simulations, the anisotropy field is taken 
Hk = 200 Oe and the saturated magnetization is M s = 
800.0 emu/cc. 



III. EXCITATION MODES IN CPP NANO- 
STRUCTURES WITH DOMAIN WALLS 

In this section, we introduce the geometry of the DW 
and study its excitation modes compared to those gen- 
erated in the uniform case. We show that the in-plane 



components of the magnetizations have distinctly differ- 
ent lowest mode frequencies that are directly related to 
the inhomogeneities of the magnetization due to the DW. 
We also discuss the magnetization evolution if we remove 
the pinning boundary conditions and allow DW to relax 
to the uniform magnetization state. In addition to the 
zero temperature dynamics, we investigate the effect of 
thermal fluctuations on the motion of transverse DW. 



Modes in the case of constrained DW 





1' 









^ M 

















FIG. 1: (online color) The CPP-SV magnetic nano-structures 
consisting of a pinned layer (PL) with fixed magnetization M p 
and free layer (FL) magnetization Mf shown for cases of ( a 
) uniform magnetization and ( b) FL with tail-to tail DW. In 
both cases, the magnetization in the middle of the free layer 
is perpendicular to the magnetization in the pinned layer. 

Figure^ shows the geometry of the systems we have 
studied. As can be seen from the schematic illustration 
in Fig. ^that in contrast with the traditional CPP-SV, 
we investigate free layer (FL) with inhomogeneous mag- 
netization due to DW in FL coupled to PL with uniform 
magnetization. We will add later in section III a third 
magnetic layer when we discuss the effect of spin torques 
on the DW. The film has an in-plane easy x-axis along 
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FIG. 2: The magnetization distribution for the three com- 
ponents of M in the DW (100 x 20 nm 2 ) as a function of 
coordinate x along the easy x-axis. 



the direction of the magnetization of the bottom pinned 
layer. The magnetization of the PL is taken homoge- 
neous. Fig. |21 shows the magnetization profile for the 
case of a small current density with DW formed in the 
plane of FL due to the uniform pinning at the x=± L/2 
boundaries, where L is the length of the side along the 
easy axis. 

Figure shows the power spectral density (PSD) 
| Mi (u) | 2 of the FL magnetization which has a peak at 
around 9.0 GHz. Before calculating the PSD, we first 
average M(x, y) over space. The size of the FL film is 
taken to be 100 x 20 x 2 nm 3 while that of the FL is 
100 x 20 x 9 nm 3 . For the unpinned boundary case, 
the average magnetization points along the easy x-axis 
and the transverse components are oscillating with a 
frequency approximately twice the ferro-magnetic res- 
onance (FMR) frequency of an infinite thin film given 

by the Kittel formula, lu = j(H k (H k + 4ttM s )) 1/2 , i.e., 
around 4.1 GHz. Figure shows that a film-size of 
1000 x 1000 x 3 nm 2 has practically the same FMR peak 
as that of an infinite thin plate. Thus, the presence of 
boundaries is an important factor which will be discussed 
throughout the rest of the paper. 

First we start discussing the DW relaxation as we re- 
move boundary pinning. This relaxation process reflects 
intrinsic DW modes. Figure|3]shows the relaxation of the 
DW when the pinning at the edges is removed at t = ns. 
The average x and z components stay zero for more than 
0.1 ns after turning off the pinning. Afterward, the x- 
component converges to M s and the z-component begins 
oscillating around zero. The y-component starts oscil- 
lating immediately around a non-zero average after the 
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FIG. 3: Spectral densities of the uniform magnetization (ar- 
bitrary units) in (a) 100 x 20 nm 2 and (b) 1000 x 100 nm 2 thin 
films as a function of the frequency (no DW) . The magnetiza- 
tion is M s = 800 emu/cc and the anisotropy is Hk = 200 Oe 
along the easy x-axis. The fluctuations of the M x component 
are small and do not show up on the scale adopted here. 



removal of the pinning at the edges and after 0.1 ns starts 
oscillating instead around zero. The initial phase of this 
decay of the DW to the uniform state shows interest- 
ing features. The x-component shows a compression- 
decompression mode which represents oscillations of the 
DW around the center and along the easy x-axis (see 
Fig. 0J>c). Simultaneously, the y-component shows a 
behavior similar to a breathing mode. Finally the last 
plot shows how a uniform magnetization which is ini- 
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tially along the hard axis relaxes to the state along the 
easy axis. Figure El shows that large damping makes the 
DW more stable to external perturbations. 



FIG. 4: Transient dynamics from the DW state to the uniform 
magnetization state with thermal fluctuations included: (a) 
magnetization components a a function of time when pinning 
has been turned off at t — ns. (b) shows the profile of the M x 
component at different time steps 1, 16, 33, 50 and 60 in units 
of At = 0.05 ns. (c) same as in (b)but for the M y component, 
(d) Relaxation of the uniform magnetization from the hard 
axis position to the easy axis position compared to that of 
the case of the DW. The damping constant is a = 0.02. The 
width of the curves in (b) and (c) represent variations in the y 
direction at a given point along the x-axis.(this figure can be 
obtained from one of the authors directly:arebei@mailaps.org 
or check the published version in Phys. Rev. B June issue 
2006) 

Magnetization dynamics of the DW can be character- 
ized in terms of its normal modes. In Fig. we show 
the spectral densities of the different components of the 
magnetization found in the DW configuration. In this 
case, the x component has a peak at lower frequency than 
the y-component. This higher frequency is directly due 
to the boundary conditions on the magnetization. The 
dependency of the magnetization on the y coordinate ap- 
pears to be very weak. As a function of the position x, 
the x and z components have a configuration which is 
odd under reflection with respect to the center, while the 
y component configuration is even. 

These different excitations of the DW can be 
qualitatively understood in a 1-D calculation with 
a simple approximation for the demagnetization field 
that of an infinite thin film. If we take, m = 
(sin#cos<p, sin#sin</?, cos 9), then the equations of mo- 
tion for the angular variables are given by 



cie 
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de 

'-—r- 
dt 

(5) 



We are looking for excitations around the ground state. 
If we take the magnetization to be in-plane, i.e, #o = ■§, 
and depends only on the x-coordinate, we find that the 
static solution should satisfy the Sine-Gordon equation 



d 2 lfn(x) 1 . . , 

j x 2 +^sin[2^ (z)] =0, 



(6) 



with the boundary condition ipo (—L/2) = n at the left 
edge and (L/2) = at the right edge. A = ^2A/K 
is the width of the DW. It should be noted here that the 
same condition arise for DW in infinite films and the only 
difference is in the boundary conditions. An analytical 
solution for this equation does not appear to be possible 
but it can be found numerically. 25 After linearization, 



e^e a + e, <p- 

0] and El become 



ipa + (/?, the equations of motion, Eqs. 



de 
It 

dp 
~dt 



2A d 2 ip(x) 2K „ dw 



2A d 2 9{x) 



Ml dx 2 
where the function G is 
2A 



G(x)t 



de 



(8) 



„, . ZA /7TV t IK 2 /nx\ 

To obtain the function G, we have Fourier transformed 
tpo( x ) an d kept only the first term. This is sufficient to 
understand qualitatively the main results of the simula- 
tion. 

The time-dependent variables e and tp satisfy ho- 
mogeneous boundary conditions and represent fluctua- 
tions around the equilibrium solution. If tp(x, t) = 
(p(x)exp(iu)t), then ip(x) can be written in the following 
form to satisfy the boundary conditions 



ip(x) = 



E 

n=l,3,5,... 



a n cos 



/ wkx \ v / mmx \ 



m=2,4,6,. 



(10) 

where x varies in the range —L/2 < x < L/2. The 
equations of motion then become algebraic equations in 
a n and b n . Then within a linear approximation, the 
magnetization components are given by 



m x 
m „ 



cos tpo(x) — (f(x) sin <po(x), (11) 
smipo(x) + <pcostpo(x). (12) 



Since the normal frequencies of the system depend on the 
wavenumber mr/L, we see that because of the parity of 
the ground state {mP x — costpo,"^ = sin ipo), the lowest 
wavenumber that appears in the y-component is larger 
(to = 2) than that of the x-component (n = 1). This 
is the reason why the breathing mode has a higher fre- 
quency than the spring (or Doring-like mode) mode, Fig. 
16171 In an infinite plane, the spring mode becomes the 
Doring mode in our case. The Doring mode is associated 
with translation of the DW, i.e., with a zero frequency 
mode. In a constrained DW, the pinned edges provide 
a restoring force and hence the center of the DW will 
oscillate instead of translating. 

The breathing mode is different from both the Doring 
and Winter modes. The Winter modes exist only in in 
infinite DW as opposed to the constrained DW treated 
here. As it can be seen from Fig. and Fig. |7|d, the 
breathing mode is the mode that is mostly excited when 
the pinned boundary conditions are turned off. 



B. Low frequency noise due to a drifting domain 
wall 

Finally in this section, we discuss the advantages be- 
hind constraining the DW to regions comparable in size 
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FIG. 5: Relaxation of the DW to the uniform state as a func- 
tion of damping for a — 0.02 and a = 0.2. 
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FIG. 6: Spectral densities of the magnetization in the con- 
strained DW (arbitrary units). The x component has a peak 
around, / = 1.0GHz. This frequency is that of the Doring- 
like mode in a finite geometry. The peak of the y-component 
is at / = 14 GHz which corresponds to the breathing mode. 
Both amplitudes for y and z have been multiplied by 100 for 
better visualization. 



to the DW width. We show that DW motion in large 
films can be the origin of l//-type noise. This noise has 
already been suggested from the measurements in Refs. 
and is detrimental to any sensing device. 
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FIG. 7: Lowest modes of the DW: The 'spring' mode (a) 
and the 'breathing' mode (b) of the DW. The dashed line 
is the equilibrium solution. The solid curves represent the 
amplitudes of the modes around the equilibrium state. Both 
x and M X:V have been normalized. 



Indeed we find that 1 / /-type noise increases if we in- 
crease the size of the thin film so it is much larger than 
the wall width A. Hence in this case, the center of the 
DW is allowed to drift away from the middle of the film 
in either direction along the easy axis due to thermal 
fluctuations and the demagnetization field. It has been 
suggested in many experiments that DW motion can give 
rise to l//-type noise 28 . Low frequency noise usually 
makes structures with non-uniform magnetization unde- 
sirable for use in magnetic sensors. However, the mag- 
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netization dynamics pattern changes dramatically if we 
constrain the DW. In the following, we first show that 
making the DW unconstrained does indeed lead to the 
l//-type noise in general agreement with experimental 
results^. 



FIG. 8: Magnetization Dynamics in the 320 x 20 nm 2 elon- 
gated CPP valve in the case of L = 320 nm S> A; A is the DW 
width : (a) The average x-component of the magnetization 
shows dynamics that gives rise to telegraph noise. The DW 
moves in a thin long strip due to random thermal excitations, 
(b) Spectral densities in the x-, -y and z-components of the 
magnetization in the elongated strip (arbitrary units). Sub- 
stantial low frequency telegraph type noise is apparent in the 
x-component. The amplitudes of the y and z peaks have been 
magnified 100 times, (c-d) Magnetization distribution in the 
xy-plane at two different times (c) and (d) that correspond 
to states with positive and negative average M x , respectively. 
This figure can be obtained from the first author or check out 
the published version in Phys. Rev. B. June 2006. 

To observe low frequency behavior in this system we 
need to reduce the effect of the restoring force on the 
DW. This amounts to enlarging the length of the sides 
along the easy axis of the film to be much larger than 
the width of the DW. This way the whole DW can move 
from the left to the right and back because of thermal 
fluctuations ( Fig. [SJ. Such behavior has been observed 
in many experiments2L2&. Figure |H1 shows the PSD's 
associated with the different components of the magne- 
tization. The x-component shows the l//-type behav- 
ior as expected. Figure |H1 shows a real-time trace for 
the average components of the magnetization. The x- 
component shows 'switching'-type behavior between two 
states, a signature of telegraph noise. The remaining 
two components are very stable and have much higher 
frequencies. The M x (t) behavior of the magnetization 
resembles the evolution of a two state system (telegraph 
noise). It is clear from Fig. [5] that this telegraph noise 
originates from the DW motion in a shallow double well 
potential. Finally, it is apparent from Fig. [HJthat DW 
magnetization appears to spend more time closer to the 
boundary than in the center. This behavior coupled with 
the telegraph-like noise are indicative of the double well 
potential for motion in the case of 'unconstrained' DW. 
In the case of constricted DW, we do not find this double 
well potential instead the potential has a single minimum 
and this is the reason behind the different noise features 
observed in both cases. However, it is not obvious what 
is the origin of such a double well potential and why it 
disappears in the case of a constricted DW. 

In the following we offer a simple explanation for such 
a behavior This difference between unconstrained and 
constrained DW's can be understood using the notion of 
charged DW introduced by Ncel 30 . The DW in elongated 
nano-elements belongs to this category and is character- 
ized by magnetic charges distributed in the vicinity of the 



DW center as schematically shown in Fig. E^a). More- 
over because of the finite size of the plate, there are signif- 
icant boundary charges at either end of opposite sign to 
that in the DW. Thus magnetostatic interaction associ- 
ated with these charges leads to a potential as a function 
of DW position as shown in Fig.^Jh). The exchange in- 
teraction contribution shown in Fig.^c) has very strong 
size dependence due to large exchange energy increase 
as DW approaches pinned boundaries. Thus, it becomes 
clear that the total potential for DW motion for uncon- 
strained DW has a double well feature which disappears 
as the DW gets constrained as shown in Fig. EJd). In 
our case the height of the barrier between the two wells 
is less than kT. 
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FIG. 9: Schematic illustration of a) magnetic charge dis- 
tribution and contributions to DW potential due to b) cor- 
responding magnetostatic energy c) exchange energy and d) 
total potential for DW motion along the easy axis. The DW, 
with width A^,, is shown to acquire negative charges due to 
V ■ M distributed around DW center and positive charges at 
the edges. The dashed line is used for the case of constrained 
DW (A™ £s L) and contrasted with the case of unconstrained 
DW (solid lines) with X m << L. 
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IV. EFFECT OF SPIN TORQUES AND ZEEMAN 
TERMS ON THE DW 

In this section we investigate the effect of spin currents 
and external fields on the magnetization in SV. 

First we study a uniformly magnetized SV and calcu- 
late the spectral density of the x-component of the mag- 
netization. We show that the spin torque can be a source 
of instabilities in this case. However in the DW case, we 
show that the effect of the spin torque can be used instead 
to control its motion. The CPP structure where this is 
possible is different than previously proposed structures. 
We instead add another magnetic layer to polarize the 
current in the direction perpendicular to the plane of the 
SV. 



A. Noise in a CPP spin valve with uniformly mag- 
netized layers 

Our discussion here will be closely related to the ex- 
perimental findings in reference ^3 where it was shown 
that spin transfer in a CPP device can give rise to 1/f- 
type noise. The noise range can be in the GHz regime 
and in effect makes the use of a CPP device as a GMR 
sensor unattractive. 

In the following we discuss a SV similar to the one 
treated in Ref. [l7] where the magnetization of the free 
layer is perpendicular to the pinned magnetization. We 
use a single spin picture to discuss the noise in this sys- 
tem. We show that this model can reproduce to a great 
extent the trend in the noise spectrum observed in the 
experiment in Ref. ^3- Adopting a single particle picture 
could be a rather crude approximation in this cas o 2 ?' 31 , 
but it is sufficient for our purpose to demonstrate the 
contribution of the spin torque to the noise of a CPP de- 
vice observed in Ref. ^3 Moreover, the single domain 
picture discussed here will help us in the interpretation 
of the numerical results of the more involved case of a 
DW. 

The CPP-SV with uniform magnetization is shown in 
Fig. We take the effective field to be equal to H off = 
(Hb — vH Cl 1500, —AttM z ) Oe where the x-component is 
much smaller than the y component. Therefore in this 
case the magnetization is expected to be almost perpen- 
dicular to the one of the pinned layer. The saturated 
magnetization is equal to 1500 emu/cc. The AFM field 
from the pinned layer is assumed small H c — 20 Oe and 
the constant v = 1 if the pinned magnetization is along 
+x and v = — 1 if the pinned magnetization is pointed 
in the — x direction. The spin torque term will be rep- 
resented by an 'effective' field term p = 1000 Oe which 
is equivalent to having a 10 mA fully polarized current 
flowing into the free layer. These parameters are chosen 
to be close to those used in the experiment of reference 
ITtI Using similar parameters, figure ITU1 shows the PSD 
for the magnetization in a 200 x 100 x 3 nm 3 thin film. 
This micromagnetic calculation clearly shows that mag- 



netization is almost uniform and is closely aligned with 
the 1500 Oe field along the y-axis. Hence we can use a 
macro-spin picture to calculate the noise spectra in this 
system. 
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FIG. 10: Spectral densities (arbitrary units) of the magneti- 
zation in the free layer in the absence of current. 

First, we need to determine the equilibrium position 
in the presence of the spin torque which is not always 
possible. The spin torque here is comparable to the pre- 
cession torque from the effective field. The equilibrium 
state is found by solving the simultaneous equations 

M x = r(H e -vH c ) (13) 

M y = r (h v - u^M z ) (14) 

M z = r (-A-KM z + v-^M y ^ (15) 

with the constraint M% + M 2 + M 2 = M 2 and r is a real 
number to be determined. These algebraic equations 
usually have up to four solutions and hence a stability 
analysis is needed to determine the stable solutions. This 
will be part of the PSD calculation of the x-component 
of the magnetization. Once, we have found the static 
solution(s) M , we make a linear expansion around it, 
M = Mo + m(i), where the perturbation is assumed 
to have the form m(t) — mexp(— iwt). The noise is 
calculated by calculating the susceptibility or the linear 
response of the magnetization due to an external small 
ac field h(i). This argument neglects the fact that es- 
tablishing a current across the layers is a non-equilibrium 
process and that a fluctuation-dissipation argument such 
as the one used below is not valid in general. However 
we have shown in Ref. [2^ that for a system in quasi- 
equilibrium, deviations from the equilibrium fluctuation 
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dissipation relation are significant only for frequencies far 
from the FMR frequency of the system. We assume in 
the following that the noise in our model depends only 
on the equilibrium state of the magnetization and hence 
only the noise around the FMR peak is well described by 
the method adopted here. 

To solve for the small perturbations from equilibrium, 
we need to solve the following system of equations, 



(iu>\ + A) • m = d, 



(16) 



where the coefficients of the matrix A are determined 
from the equations of motion for the magnetization, 



in the simulations^! and the experiment, we may be able 
to deduce the switching indirectly since depending on the 
value of the H e field, we may end up with more than one 
possible solution to the static equations. For large neg- 
ative easy axis fields (Fig. 1111) . we see the usual shape 
of FMR curves. The PSD curves in this section only 
are normalized differently from those in other sections of 
the paper. The damping parameter in this calculation 
is taken a = 0.005, which is appropriate for a permalloy 
even though we expect a higher value due to spin accu- 
mulation at the interfaces between a normal conductor 
and a ferromagnet. 20 



A n = 

A12 = 

A13 = 

di = 



(2H e Mo tV — H y Mo lX ) 



+ W (~ 2v V M ^y + 4ttM s m 0}Z ) , 

-jjji {H e M , z + 4ttM 0;X M , z ) 

+ W ( _2i/pMo >* + A7rM s M °,y + H v M s) (17) 
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{M 2 



2 y + Ml z ), 



The remaining coefficients can be determined in a similar 
way. 

The coefficients in front of the x-component of the 
ac field h(t) are grouped in the vector d. The stable 
solutions will be those for which the imaginary part of 
u> is negative or zero, det \ioj + A| =0. In the absence 
of the spin torque, the frequencies are real in a stable 
system. The imaginary frequencies that appear are a 
signature that the spin torque can act as a (damping) 
force. The noise spectrum is found by solving for m in 
Eq. 1161 In the experiment only the noise in x-component, 
C xx (w) = J dt{M x (t)M x (0))e iuJt along the pinned mag- 
netization is of interest. It is found from the fluctuation- 
dissipation relation at inverse temperature /3 



C xx H = icoth 



Im- 





d x 


Al2 


A13 


det 


d 2 


iuj + A 2 2 


A 23 




d 3 


A 3 2 


ILU + A 33 



det \iu) + A| 



(18) 

These steps are carried out for all the static solutions 
that are found for each bias field H e in the presence of 
the spin torque. 

The magnetization of the pinned layer is taken in the 
-x direction (as in the experiment) and the current is pos- 
itive when it flows from the pinned to the free layer. In 
this case we expect to see more noise for negative easy 
axis fields and less noise for positive easy axis fields. The 
l//-type noise is observed when the field along the easy 
axis is small and negative. Since this equilibrium analy- 
sis can not show actual switching between two states as 




FIG. 11: The FMR curve for the free layer with uniform 
magnetization and for negative bias fields (arbitrary units). 

In figure IT21 we plot the noise for H e = 100 Oe and 
H e = —100 Oe. Clearly for the case with the positive 
field, the noise is completely suppressed compared to the 
case with negative easy axis biasing. This is consistent 
with the experiment. Therefore the state with positive 
biasing is equivalent to a state with large effective damp- 
ing. This large damping is coming from the spin mo- 
mentum transfer. If we turn off the current, we get back 
the usual FMR (bright) spectrum (see Fig. I13|) in this 
case too. This asymmetry between positive and negative 
biasing fields close to the perpendicular direction of the 
free layer will be important later when we have a DW in 
the presence of spin torques. 

Therefore the single spin model captures the 'bright' 
and 'dark' regions of the spectral density for frequencies 
around the FMR frequency (see Fig. 2 in Ref. nm . Fig. 
114b shows the strength of the power as a function of the 
negative bias field. Clearly for large biasing we have 
less noise as expected. Now, if we plot the same curve 
for positive fields, Fig. 114b . we find a very interesting 
result. For large positive fields, we have the usual 'dark' 
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stable than the bright one in this case. The x-component 
of the magnetization in the dark state is negative, i.e., 
opposite to the direction of the easy axis field while the 
bright one is along the field H e . This is most probably 
the origin of the l//-type noise in the system. The 1/ / 
region in the experiment appears on the negative side 
of the easy axis field. Here it appears on the positive 
sidsit The reason is that the zero point of the axis 
is not well known in the experiment. The experiment 
estimates that the magnetization is perpendicular to the 
pinned layer at H e = 34 Oe which should correspond to 
—20 Oc in our case. Therefore there is a shift of about 
50 Oe in the reference point which is approximately the 
field when two states become possible as a solution to our 
equations. The important point we need to remember 
that the spin system behaves differently for positive and 
negative bias when there is a spin torque. This is mainly 
due to the fact that in one case the spin torque is acting 
as a regular field, while in the other, it is acting as an 
extra source of damping. 



FIG. 12: Comparison of the PSD (arbitrary units) for positive 
and negative biasing fields H e . For H e = 100 Oe the peak 
(which is multiplied by 500) has been completely suppressed 
and there is a shift to the left, a signature of an over-damped 
state. 



B. Tri-layer CPP structure with a trapped DW 

Next, we turn to the study of the DW case. First, we 
show how a DW in CPP-SV can be manipulated by low 
currents through the spin torque. The interaction of the 
DW with an external field will be also shown. 




2.0 4.0 6.0 
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FIG. 13: Effect of the spin torque on the PSD (arbitrary 
units) curve for H e = 100 Oe. The spin torque is clearly 
acting as a damping force for positive current. For 1 = 0, 
the PSD curve has the familiar FMR shape. 



regions that reflect high damping states. As we lower the 
field, we find that the system now can sustain two states, 
one bright and one is dark. The dark state is actually less 



1. The effect of spin polarized current 

First, we consider an alternative CPP structure. As 
will be shown below this modification of the traditional 
CPP structure can be done for at least three reasons. 
One reason is to create structures where DW can be 
manipulated with spin momentum in most efficient way. 
Secondly, we would like to be able to detect domain wall 
motion with GMR effect. We also find that the suggested 
CPP structure modifications may have some advantages 
in terms of reducing effects of magnetization instabilities 
due to the spin momentum transfer. As has been shown, 
for example in Ref. l3ll even in CPP devices with nom- 
inally uniform magnetization, the spin torque can give 
rise to magnetization instabilities. Thus, the latter rea- 
son should be kept in mind as an important one. 

In the following we consider a CPP structure which 
has three magnetic layers (Fig. I15fl where the DW layer 
is sandwiched between two pinned magnetic layers with 
one of them polarized along the direction of the current 
and the other polarized along the easy axis of the middle 
DW layer. In the following simulations, the bottom 
layer is taken to be 100 x 20 x 9 nm 3 , the middle layer 
is 100 x x20 x 2 nm 3 , and the reference layer has the 
dimensions 100 x 20 x 3 nm 3 . In this geometry, the two 
outer magnetic layers lead to a two different spin torques, 
r = V b + IV, acting on the middle magnetic layer 

r = —pi (m x (m x nig) — m x (m x m^)) , (19) 
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FIG. 15: (Online color) Schematic geometry layout for the 3 
magnetic -layers CPP structure. The first layer is now polar- 
ized along the current direction. The middle has the DW. The 
third layer has been introduced to enable the GMR sensing 
of the DW motion. 
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FIG. 14: PSD (arbitrary units) for both negative and posi- 
tive biasing: (a) PSD (FMR) for negative easy axis biasing 
fields. (b) PSD curves for positive bias fields. For fields 
approximately between 50 Oe and 20 Oe, there are two possi- 
ble states for the system; one is over-damped (dark) and the 
other is regular (light). Outside this range of fields, only dark 
or light states exist. 



where ( m-r) is the magnetization direction of the 
bottom (top) layer. The damping parameter in this 
section has been increased to a = 0.08 to better account 
for spin accumulation^. 

First, we investigate if the DW can be moved along 



the easy x-axis with moderate currents. This will en- 
able a spin torque Tb = — m x h^, with an effective field 
along the x-axis and proportional to the y-component of 
the magnetization in the DW layer that is largest at the 
center, h^, = pl(m y , — m x , 0). Since M y « M s around 
x = 0, this gives us the optimal field needed to push the 
DW off the center and this appears to be a primary rea- 
son why only very low currents are needed to have an 
appreciable motion of the DW in considered CPP geom- 
etry. FigureED shows the effect of the spin torque on the 
DW in the 3-layer geometry as a function of the current. 
We find that the spin torque from the top layer has a rel- 
atively small effect on the dynamics of the DW since its 
effective field is = pl(0, —m z ,m y ). Given that the 
z-component of the magnetization is practically zero for 
the currents in the case shown in Fig. ^| the effect of 
on the magnetization is negligible. We find that indeed 
in this geometry, the spin torque can be used to con- 
trol the motion of DW with very small current densities. 
This is primarily due to the fact that the constrained DW 
has a non-zero y-component of magnetization in the DW 
region. 

The displacement of the DW by the spin torque (we 
make sure that the Oersted field is not the origin of this 
motion) is easily understood from the equation of motion 
without the demagnetization field. Taking account of 
only the spin torque, the exchange and the anisotropy, 
the static equations for the magnetization are 



2A d 2 m x 2K 



cm x , 



2A d 2 m y 
MF dx 2 



pl(m x 



2 A d 2 m z 
M? dx 2 



■ plrriy 



an , 



(20) 
(21) 
(22) 



11 



where c(x) is a real function and m 2 = 1. Clearly, in 
this case the x-component is coupled to the y component 
which acts as a source term for the x-component. Ne- 
glecting anisotropy and integrating the equation for the 
m x component around zero, we find that the difference 
in the slope of m x (x) for x=0 and small x = e is given 
by 



dm 2 
dx 



x=Q 



dm x 
dx 



pleniy. 



(23) 



For positive e, i.e. a shift to the right, the slope at x = is 
smaller than that at x = e which is approximately equal 
to that at x = in the absence of current. Therefore 
/ should be negative for positive m y which is approxi- 
mately equal to the one around x = 0. This is confirmed 
by the numerical integration of LL equation in fig. 1161 
The spin torque can therefore be used to move the DW 
in a controlled fashion with low currents. In addition we 
find, that a three-layer structure may actually have lower 
frequency noise in the presence of spin torque than the 
structure investigated in Ref. I'ilL This potential advan- 
tage of our proposed structure is however realized only if 
the middle layer geometrical dimensions are comparable 
with the DW width, Fig. [T7J In this case the PSD in 
the x-component does not show any substantial low fre- 
quency noise. For the parameters used here, we find that 
the DW width is approximately 40 nm. The dimension 
of the film is 100 nm. Therefore the DW is barely con- 
strained and hence the reason behind the sensitivity of 
the DW to external forces due to fields or currents. 
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FIG. 16: Magnetization profiles of the in-plane components as 
a function of the sign of the current in the CPP device. In (a) , 
depending on the sign of the current the curve shifts to the 
left or to the right, b) The shape of the x and y components 
of the magnetization in the DW is a slight distortion from 
that without currents. 
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FIG. 17: PSD (arbitrary units) of the different components 
of the magnetization in the presence of current / = 0.001 mA 
for a constrained DW. 
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FIG. 18: The time evolution ol the different components of 
the average magnetization of the DW for / = 0.1 mA. (a,c) 
and I = —0.1 mA (b,d). No stationary solutions exist at high 
currents. 
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At higher currents, we are no longer in a linear regime. 
The numerical integration of the LLG equation shows 
that the z-component of magnetization becomes more 
significant as we increase the current and this contributes 
to the twist of the DW and no stationary solutions are 
possible in this case. FigurelT%lshows the time evolution 
of the magnetization for current densities of the order of 
1.0 x 10 7 A/cm 2 . 

The magnetization dynamics is a regular periodic ro- 
tation. In this case the spin torque can be used to selec- 
tively excite higher modes of the magnetization as com- 
pared to those studied in section II. 

2. The Effect of External Magnetic Field on a DW 

Finally in this section, we investigate effect of an ex- 
ternal magnetic field. We add a Zeeman term to the total 
energy and study the displacement of the DW due to an 
external field along the easy axis. 

The external field along the easy axis is applied to the 
middle layer in the presence of a small current to mea- 
sure resistance changes across the CPP structure and so 
that no spin torque effects are appreciable on the DW. 
Interestingly the three-layer structure with DW does not 
require biasing which is needed for standard CPP struc- 
tures to achieve linear dependence of resistance on the 
external field. The calculated transfer curve of resis- 
tance R versus field H is shown in Fig. ^] As can 
be seen this dependence is centered around zero, has a 
large slope dR/dH, and shows small hysteresis. More- 
over, the system appears to be more stable to pertur- 
bations by the spin torque and no 1/f-type behavior is 
observed in this case. Our device is therefore well suited 
to function as a magnetic sensor. However the proposed 
structure lacks an important property which is needed 
in memory applications and that is non-volatility. As 
we remove the voltage across the CPP-SV, the DW re- 
laxes back to its equilibrium position and hence any state 
stored in the DW position is lost. Nevertheless, our pro- 
posed CPP-SV structure can be incorporated as part 
of a logic device. Recently, properly redesigned CPP 
structures have been proposed for reprogammable logic 
elements22ii&. This latter application does not require 
non- volatility and henc e our device can be utilized in a 
similar way as in Refs. I.Sfl.Sfl Our device has an advan- 
tage compared to that proposed in Ref. and that is 
only much smaller currents are needed in our case. 
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FIG. 19: (a) Resistance versus external field applied along 
the easy x axis for the 100 x 20 nm 2 film. The reference layer 
is pinned in the +x direction. The resistance is normalized 
in such a way that for R = 1 the two layers are parallel 
and for R = 1.1 they are anti-parallel, (b-d) Profiles of the 
magnetization in the xy plane for zero (b) and positive fields 
(c-d), respectively. The scale for the x and y axis are different. 
For a given x, five points are plotted along the y-axis. The 
DW is fairly uniform along the y axis. 
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V. SUMMARY AND DISCUSSION 

In summary, we have presented a study of magnetiza- 
tion dynamics for CPP geometry which includes a con- 
strained DW layer. We have identified a Doring-type 
mode and a new breathing mode. It is shown that the 
lowest modes of the DW dynamics can be understood in 
terms of the parity of the inhomogeneous ground state. 
We have investigated in details how the constrained DW 
dynamics is affected by the spin polarized current and 
thermal fluctuations and compared it with the traditional 
single domain free layer structures. In particular, we 
found that the currents needed to measure any apprecia- 
ble motion of the DW are at least two orders of magni- 
tude less than usual values of currents needed to switch 
the single domain magnetization. This difference is at- 
tributed to the appearance of a significant magnetization 
component of the constrained DW that is perpendicular 
to the pinned layer magnetization and the exchange field. 

We also find that thermally activated motion of the 
constricted DW has lower weight in the lower frequency 
region than that of the unconstrained DW. The latter 
shows well known telegraph type noise characteristics. 
This difference can be understood using notion of charged 
DW introduced by Neel^S and competition of the mag- 
netostatic and size dependent exchange interaction con- 
tributions to the DW potential for motion along the easy 



axis (see Fig. |5J) . The three magnetic layer CPP struc- 
ture was introduced so that the spin torque effect on DW 
layer is maximized . This CPP geometry has been inves- 
tigated and found to have a number of interesting prop- 
erties such as (i) DW can be easily controlled by an exter- 
nal field or a polarized current with relatively small cur- 
rent densities; (ii) linear dependence of resistance on the 
external field and current; (iii) improved magnetization 
stability characteristics. Experimental realization of the 
device proposed here requires finding ways to constrain 
the DW within the middle 'free' layer. The pinning at 
the edges can be realized by creating permanent magnets 
with different coercivity and/or with anti- ferromagnetic 
(AF) coupling. The AF coupling at the edges needs anti- 
ferromagnets with different Neel temperatures^ so that 
a properly designed field-cooling procedure could lead to 
the pinning in opposite directions at the edges of the 
magnetic stripe. Other alternatives such as special shap- 
ing and padding also have been discussed in the litera- 
ture in the context of stability of DW in magnetic nano- 
elements^. 
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